Universal Fluctuations of Growing Interfaces: Evidence in Turbulent Liquid Crystals 
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We investigate growing interfaces of topological-defect turbulence in the electroconvection of ne- 
matic liquid crystals. The interfaces exhibit self-afiine roughening characterized by both spatial and 
temporal scaling laws of the Kardar-Parisi-Zhang theory in 1+1 dimensions. Moreover, we reveal 
that the distribution and the two-point correlation of the interface fluctuations are universal ones 
governed by the largest eigenvalue of random matrices. This provides quantitative experimental 
evidence of the universality prescribing detailed information of scale-invariant fluctuations. 
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Growth phenomena have been a subject of extensive 
studies in physics and beyond, because of the ubiquity 
in nature and of their importance in both engineering 
and fundamental science. Over decades, physicists have 
found that growth phenomena due to local processes typ- 
ically lead to the formation of rough self-affine interfaces, 
as exemplified in paper wetting, burning fronts, bacterial 
colonies, and material morphology, to name but a few, 
and also in various numerical models jl| . Being obviously 
irreversible, local growth processes provide a challenging 
situation toward understanding the scale invariance and 
the consequent universality out of equilibrium. 

The roughness of interfaces is often quantified by their 
width w{l,t) defined as the standard deviation of the 
interface height h{x, t) over a length scale I at time t. 
The self-affinity of interfaces then implies the following 
Family- Vicsek scaling 
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with two characteristic exponents a and /3, the dynamic 
exponent z = a//?, and a crossover length scale Z* t^l^ . 

The simplest theory to describe such local growth pro- 
cesses was proposed by Kardar, Parisi, and Zhang (KPZ) 
on the basis of the coarse-grained stochastic equation 
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dt 



h{x, t) = vo + vV^h + -{Vh f + £,{x, t) (2) 



with {£,{x,t)) = and {£,{x,t)i{x' ,t')) = D5{x - x')6{t - 
t'). For l-f-l dimensions, the renormalization group ap- 
proach provides exact values of the exponents at a^^^ = 
1/2 and /3^P^ = 1/3 which are universal as widely 

confirmed in numerical models [ij. Moreover, the (1-1-1)- 
dimensional KPZ class attracts growing interest thanks 
to rigorous work on the asymptotic form of the fiuctu- 
ations in solvable models |^-l6|. This opens up a new 
aspect in the study of scale-invariant phenomena toward 
the universality beyond the scaling laws. 

In contrast with such remarkable progress in theory, 
the situation in experiments has been quite different. A 
considerable number of experiments have been performed 



on various growth processes f[\ and confirmed the ubiq- 
uity of rough interfaces. Concerning the universality, 
however, experimentally measured values of the expo- 
nents are widely diverse and mostly far from the KPZ 
values for both a and /3 [ll|. To our knowledge, only two 
experiments among dozens directly found the KPZ ex- 
ponents: in colonies of mutant bacteria [3] and in slow 
combustion of paper [sj, apart from few other experi- 
ments showing indirect indications [^, [lo| . One of the 
main difficulties shared by most experiments, including 
the above two, is that one needs to repeat a large num- 
ber of experiments in the same controlled conditions to 
accumulate sufficient statistics. In this Letter, studying 
growing interfaces of turbulent liquid crystals, we over- 
come this difficulty and report clear experimental evi- 
dence of not only the universal scaling laws but also the 
universal fluctuations of the KPZ class through critical 
comparisons with the wealth of theoretical predictions. 

The electroconvection occurs when an external volt- 
age is applied to a thin layer of nematic liquid crystal, 
triggering the Carr-Helfrich instability 



11 



We focus 

on interfaces between two topologically different turbu- 
lent states called the dynamic scattering modes 1 and 2 
(DSMl and DSM2), which are observed with sufficiently 
large voltages. The essential difference between them 
lies in the density of topological defects called the discli- 
nations. Upon applying a voltage, we first observe the 
DSMl state with practically no defects in the director 
field, which lasts until a disclination is finally created 
owing to the breakdown of surface anchoring This 
forms a DSM2 cluster composed of a large quantity of 
disclinations, which are constantly elongated, split, and 
transported by fiuctuating turbulent flow around. While 
DSM2 may coexist with DSMl in a regime of spatiotem- 
poral intermittency [l5], for larger voltages we observe 
growing DSM2 clusters driven by the above-mentioned 
stochastic local contamination processes. 

Our experimental setup consists of a quasi-two- 
dimensional sample cell, an optical microscope, a ther- 
mocontroUer, and an ultraviolet pulse laser (see Rcf. 
[l3| for detailed descriptions). The cell is made of 
two parallel glass plates with transparent electrodes. 
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FIG. 1: (Color online) Growing DSM2 cluster, (a) Images. 
Indicated below is the elapsed time after the emission of laser 
pulses, (b) Snapshots of the interfaces taken every 5 s in the 
range 2s < t < 27s. The gray dashed circle shows the mean 
radius of all the droplets a,t t — 27 s. The coordinate x at this 
time is defined along this circle. 



which are spaced by a polyester film of thickness 12 /im 
enclosing a region of 16 mm x 16 mm for the convec- 
tion. We chose here the homeotropic alignment of liquid 
crystals in order to work with isotropic DSM2 growth, 
which is realized by coating Af,iV-dimethyl-iV-octadecyl- 
3-aminopropyltrimethoxysilyl chloride uniformly on the 
electrodes using a spin coater. The cell is then filled with 
A'^-(4-methoxybenzylidene)-4-butylaniline doped with 
0.01 wt.% of tetra-n-butylammonium bromide. The cut- 
off frequency of the conductive regime ll| is 850 ± 50 Hz. 
The cell is maintained at a constant temperature 25.0 °C 
with typical fluctuations in the order of 10"'^ K. The con- 
vection is observed through the transmitted light from 
light-emitting diodes and recorded by a CCD camera. 

For each run we apply a voltage of 26 V at 250 Hz, 
which is sufficiently larger than the DSM1-DSM2 thresh- 
old at 20.7 V. After waiting a few seconds, we shoot into 
the cell two successive laser pulses of wavelength 355 nm 



and energy 6nJ to trigger a DSM2 nucleus [13j. Figure 
[T] displays typical growth of a DSM2 cluster. We repeat 
it 563 times to characterize the growth process precisely. 

We define the local radius R{x, t) along the circle which 
denotes the statistically averaged shape of the droplets, 
as sketched in Fig. [IJb). This measures the interfacial 
width w{l,t) = (^y{[Rlx,t) - {R)i]^)i) and the height- 
difference correlation function C{l,t) = {[R{x + l,t) — 
R{x, t)]'^), where (• • and (• • •} denote the average over 
a segment of length I and all over the interface and ensem- 
bles, respectively. Both w{l, t) and C{1, tY^'^ are common 
quantities for characterizing the roughness, for which the 
Family- Vicsek scaling [Eq. dT])] is expected. 

This is tested in Fig. [51 Raw data of w{l,t) and 
C{l,ty^^ measured at different times [Fig. [2l^a,b)] grow 
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FIG. 2: (Color online) Scaling of the width w{l,t) and the 
height-difference correlation function C(l,t). (a,b) Raw data 
of w{l, t) (a) and C{1, tf^ (b) at different times t. The length 
scale I is varied up to 2tt{R) and it{R}, respectively, (c) 
Time evolution of the overall width W{t) and the plateau 
level Cpi(t)^^^ of the correlation function, (d) Collapse of the 
data in (a) showing the Family- Vicsek scaling [Eq. ([T])]. The 
dashed lines are guides for the eyes showing the KPZ scaling. 



algebraically for short length scales I ^ l^, and converge 
to constants for I ^ l^, in agreement with Eq. ([T]). The 
power a of the algebraic regime measured in the last 
frame t — 28.4s is found to be a = 0.50(5). Here, the 
number in the parentheses indicates the range of error 
in the last digit, which is estimated both from the un- 
certainty in a single fit and from the dependence on the 
fitting range. The found value of a is in good agreement 
with the KPZ roughness exponent a^^^ = 1/2. 

The temporal growth of the roughness is measured by 
the overall width W{t) = ^J{[R{x,t) - (R)]^) and the 
plateau level of the correlation function, Cpi(t)^/^, de- 
fined as the mean value of C{l,ty^^ in the plateau re- 
gion of Fig. [2][b). Both quantities show a very clear 
power law t>^ with = 0.336(11) [Fig.djc)] in remarkable 
agreement with the KPZ growth exponent = 1/3. 

Furthermore, rescaling both axes in Fig. [IJa) with the 
KPZ exponents, we confirm that our data of w{l,t) col- 
lapse reasonably well onto a single curve [Fig. [2jd)]. A 
collapse of the same quality is obtained for C{l,ty^^. 
We therefore safely conclude that the DSM2 interfacial 
growth belongs to the (l-l-l)-dimensional KPZ class. In 
passing, this rules out the logarithmic temporal scaling 
claimed by Escudero for the droplet geometry 14 1. 

Our statistically clean data motivate us to test further 
predictions on the KPZ class beyond those for the scaling. 
In this respect one of the most challenging benchmarks 
may be the asymptotic distribution of height fluctua- 
tions, calculated exactly for solvable models A gen- 
eral expression was proposed by Prahofer and Spohn \^ , 
which reads h{t) ~ Vaot + iA^Xt/2y/^x with A = D/2v, 
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FIG. 3: (Color online) Parameter estimation, (a) Growth rate 
d{R)/dt averaged over 1.0 s against t~^^^ . The {/-intercept of 
the linear regression (dashed line) provides an estimate of A. 
(b) C{l,t)/l against I for different times t. Inset: nominal 
estimates of A obtained from w{l, t) (blue bottom symbols) 
and C(l, t) (green top symbols) as functions of t (see text). 



the asymptotic growth rate Voa, and a random variable 
X obeying the Tracy- Widom (TW) distribution fl5|, or 
the (rescaled) largest eigenvalue distribution of large ran- 
dom matrices. The random matrices are from the Gaus- 
sian unitary and orthogonal ensemble (GUE and GOE) 
(l6j for curved and flat interfaces, respectively. This im- 
plies an intriguing relation to the random matrix the- 
ory and requires no fitting parameter provided that the 
values of the two KPZ parameters A and A are mea- 
sured. The prediction was tested once for flat interfaces 
in the paper combustion experiment [l7| with an appar- 
ent agreement. However, the authors had to shift and 
rescale the distribution function for want of the values of 
the KPZ parameters, in which case the difference among 
the predicted distributions and the Gaussian one is un- 
pronounced. They also had to discard data subject to 
intermittent advance of burning fronts due to quenched 
disorder Therefore, a quantitative test of Prahofer 
and Spohn's prediction has not been carried out so far. 

We first measure the value of A experimentally. For 
the circular interfaces, A is given as the asymptotic radial 
growth rate, which has a lea ding correction term as A ~ 
— > oo Il8l . This relation is indeed 
confirmed in Fig. [3l^a) and yields a precise estimate at 
A = 35.40(23) Aim/s. 

The parameter A can be determined, at least for fiat in- 
terfaces, from the amplitude of C(/, t) and w{l, t) through 
C ~ a; and w"^ ~ Al/Q in the limit t oo 18]. Fig- 
ure [S^b) shows C{l,t)/l against / for different times t. A 
similar series of plots is obtained for 6w^/l. The value 
of A can be estimated from the plateau level or the lo- 
cal maximum of these plots, but we find that these es- 
timates increase slowly with time and do not agree with 
each other (inset). This allows us to have only a rough 
estimate A w 10 /im for the range of time we study. 

Now we test Prahofer and Spohn's prediction for the 
circular interfaces: 
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FIG. 4: (Color online) Local radius distributions, (a) Cumu- 
lants {R")c vs t. The dashed lines are guides for the eyes 
showing the indicated powers, (b) Skewness {R^}c/ {R^fJ^ 
and kurtosis {R"^) c / {R^)1 ■ The dashed and dotted lines indi- 
cate the values of the skewness and the kurtosis of the GUE 
and GOE TW distributions, (c) Local radius distributions 
as functions oi q = {R - \t) / {A^\t/2f^ . The dashed and 
dotted lines show the GUE and GOE TW distributions, re- 
spectively, (d) Differences in the cumulants of q and xgue- 
The dashed line indicates (g")c = (Xgue)c- Inset: the same 
data for n = 1 in logarithmic scales. The dashed line is a 
guide for the eyes. 



with a random variable xgue obeying the GUE TW dis- 
tribution. We first compute the cumulant (i?")c, for 
which Eq. dSl) implies (i?")c ~ {A^\/2Y/^ {xqv^)^^^'^ 
for 71 > 2. Our data indeed show this power- law be- 
havior in time [Fig. lUJa)], though higher order cumu- 
lants are statistically more demanding and hence provide 
less conclusive results. We then calculate the skewness 
{R^)J{R^)V^ and the kurtosis {R^)J{RYc, which do 
not depend on the parameter estimates. The result in 
Fig. m^b) shows that both amplitude ratios asymptoti- 
cally converge to the values of the GUE TW distribution, 
about 0.2241 for the skewness and 0.09345 for the kurto- 
sis 0], and clearly rules out the GOE TW and Gaussian 
distributions. Conversely, if we admit the GUE TW dis- 
tribution, the amplitude of {R^)c offers a precise estimate 
of A at 9.98(7) //m, which is consistent with the direct es- 
timate obtained above and hence used in the following. 

Histograms of the local radius R{x,t) are then made 
and shown in Fig. Hl[c) for two different times as func- 
tions oi q = {R - \t)/[A^\t/'2.)^/^, which corresponds 
to Xgue if Eq. ([3]) holds. The experimental distributions 
show remarkable agreement with the GUE TW one with- 
out any fitting, apart from a slight horizontal translation. 
Indeed, time series of the difference between the nth or- 
der cumulants of q and xgue [Fig. IDJd)] reveal that the 
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FIG. 5: (Color online) Two-point correlation function C2{l,t) 
plotted in the rescaled axes u = {Al/2){A^ Xt /2)''^^^ and 
C2(u) = C2{l)/{A^\t/2f/'-^. The dashed line indicates the 
theoretical prediction g2(u). Inset: integral of the rescaled 
correlation function = C2{u)du as a function of time 
t. The difference from the theoretical counterpart (72"* = 
/o°° g2(u)Au is shown. The dashed line is a guide for the eyes. 

second to fourth order cumulants of q converge quickly 
to the GUE TW values, while the first order one, i.e., the 
mean, algebraically approaches it with a power close to 
— 1/3 (inset). This is theoretically reasonable behavior 
which stems from the existence of an additional constant 
term in Eq. ([3]). Therefore, we conclude that the local 
radii of the DSM2 nuclei asymptotically obey the GUE 
TW distribution at least up to the fourth order cumu- 
lants, confirming Prahofer and Spohn's prediction. 

We also measure the two-point correlation function 
C2{l,t) = {R{x + l,t)R{x,t)) - (R)^. Theory predicts 
that 02(1, t) is asymptotically described by the Airy2 
process A2{t) or by the dynamics of the largest eigen- 
value in Dyson's Brownian motion of GUE matrices [l6| 
as C2{l,t) ~ {A^Xt/2)^/^g2iu) with g2(u) = {A2{u + 
t)A2{t)) and u = (AV2)(A2At/2)-2/3 0. Our exper- 
imental data confirm this with an algebraic finite-time 
correction consistent with the power —1/3 (Fig.[5|). 

In comparison with past experimental studies showing 
diverse scalings, one may wonder why the liquid crys- 
tal turbulence exhibits such clear KPZ-class behavior. 
We consider that the following three factors are essen- 
tial: (a) The growth of DSM2 results from strictly local 
processes due to the turbulent flow on the interfaces and 
not from inward or outward interactions of the cluster, 
which could induce long-range effects and affect the uni- 
versality, (b) The stochasticity of the process stems from 
intrinsic turbulent fluctuations overwhelming quenched 
disorder, (c) Good controllability and fast response of 
the liquid crystals allowed us to repeat hundreds of ex- 
periments in the same conditions, leading to statistically 
reliable data. The reproducibility of the presented re- 
sults was confirmed with different voltages and spatial 
resolutions with the same quality of data (not shown) . 

In conclusion, measuring the growth of DSM2 nuclei in 
the electroconvection, we have found the circular inter- 
face roughening clearly characterized by the scaling laws 
of the KPZ class in 1-1-1 dimensions. Moreover, we have 



shown without fitting that the fluctuations of the cluster 
local radius asymptotically obey the Tracy- Widom dis- 
tribution of the GUE random matrices and revealed the 
finite-time effect. Together with the agreement in the 
two-point correlation, our experimental results quantita- 
tively confirm the geometry-dependent universality of the 
(l-l-l)-dimensional KPZ class prescribing detailed infor- 
mation of the scale- invariant fluctuations. In this respect, 
investigations of flat interfaces in the same system are of 
outstanding importance and are in progress. 
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us with numerical values of the TW distributions and 
the covariance of the Airy2 process. This work is partly 
supported by JSPS and by MEXT (No. 18068005). 

Note added in proof. After submission of this Let- 
ter, Sasamoto and Spohn reported an exact solution of 
the (l-l-l)-dimensional KPZ equation fio*], which offers a 
clear theoretical ground of our experimental results. 
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